Inaccurate use of asymptotic formulas 
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The asymptotic form of the plane wave decomposition into spherical waves, which is used to 
express the scattering amplitude in terms of phase shifts, is incorrect. We explain why and show 
how to circumvent the mathematical inconsistency. 

In quantum mechanics the foUowing plane wave decomposition into spherical waves is used 
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where k ■ r = fcr cos©, Pi is the ^th Legendre polynomial, and ji is the Ith spherical Bessel function. Because we are 
interested in the large distance behavior of the wave function in scattering theory, we need the asymptotic form of 
the spherical function 



and we rewrite Eq. ^ as 
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3'^ '' I Vi'(2Z + l)P,(cose)sin(fcr-7rV2). (3) 



Equation iPJ is given in numerous textbooks on quantum mechanics, including Schiff^ and Landau and Lifshitzi^ 
Astonishingly, the expression ^ is meaningless, and for this reason we put the question mark over the approximate 
equality. The series is not only divergent, but it cannot even be treated as an asymptotic expansion of the function 
e^^ '' at large distances. 

To see how badly the series Q diverges, we consider the special case cos 9 = 1. Then, — 1, and after a simple 

calculation we obtain 



sin(fcr) .cos(fcr) 



y(2Z + l)cos2(7rV2)-i— -^y(2Z + l)sin2(7rZ/2) (4a) 
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kr E(4" + l)-=^E(4« + 3). (4b) 
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According to Eq. 10}, both the real and imaginary parts of e**''' contain a divergent series, for any value of r. 

However, we often consider asymptotic series that are divergent, but still correctly represent certain functions. The 
infinite series ao(x) -\-ai{x) +a2(x) + . . . is the asymptotic expansion of the function /(x) at xq (which can be infinite) 
i^ 
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~rT(/(^) ~ E"'(^)) ^ ^ for X ^ Sq. (5) 
Equivalently, the series is asymptotic if 



Unix) 
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> (J tor a; ^ Xq. (d) 

ai(x) 

Due to the definition |SJl, any finite subseries of an asymptotic series approximates the function f{x) and the ap- 
proximation becomes better and better as a; — > Xq. However, the series ^ does not satisfy the condition (jSJ, and 
consequently, it cannot be treated as an asymptotic expansion of e*"^ "" at large distances. 
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What is wrong with the expansion It appears that the approximate formula ^ requires that 

fcr»i;G + l). (7) 

For completeness, we derive this condition here, and find not only the first but also the second term of the 1/z 
expansion of ji{z). It is well known (see, for example, Ref. 0) that the spherical Bessel functions can be written as 



, , , , 1 d \ ' sin z 



If we use Eq. © and the recursion formula, 

= (9) 

we can easily prove by induction that 

, , sinfz — 7r//2) 1 , , cos(z — 7r//2) , 1 , , , 

Mz) = \ ^ ' + + 1) ,2 + (10) 

If we compare the two terms of the expansion pU|l . we find that the approximation ||2Jl holds if the condition {Tj) 
is satisfied. When we perform the summation in Eq. JSjl, we find that the terms for sufficiently large I violate the 
requirement and effectively destroy even the approximate equality. 

Although the decomposition Q is incorrect, the results obtained by means of it are usually correct. Obviously, the 
famous formula, that expresses the scattering amplitude via the phase shifts, is correct. However, it is of interest to 
see why the derivation works. Therefore, we first discuss the standard procedure, which can be found, for example, 
m Refs.QandU and then we show how to avoid the mathematical inconsistency. 

By assuming azimuthal symmetry, the scattered wave function is 
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0k(r) =^A,i'(2/ + l)P,(cose)i?i(r), (11) 
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where the Ri are the radial wave functions and the Ai are coefficients to be determined. If we assume that the 
asymptotics of the radial functions are 

kr 

where 6i denotes the /th phase-shift, we can rewrite Eq. as 

? 1 

0k(r) « — VA,i'(2; + l)P,(cose)sin(fcr-7r//2 + (5/). (13) 
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We still put the question mark over the equalities that are mathematically inappropriate. 

Now, we compare the wave function H13|) with the expected asymptotic form of the scattered wave function 

0k(r)=e'''"- + /(e)^, (14) 
where /(B) is the scattering ampHtude. If we use the plane-wave decomposition JSJ, we find the equation 
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« i'(2/ + l)Pi(cose)sin(fcr-7r//2) + /(e) — . (15) 
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If we equate the terms proportional to e , we find that Ai = e'*' , which when substituted into the terms proportional 
to e*'^'", provides the well known result 

/(Q) = ^E(2^ + l)^Kcose)[e2'^' - 1]. (16) 

* (=0 
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We obtained the relation H16() using the mathematically meaningless equations Q, H13() . and H15() . Next, we show 
how to derive Eq. (|16|l avoiding the inconsistency. 

Again we start with the wave function in the form 111(1 . and we use the asymptotics of the radial wave function 
(|12|l . but only for fixed values of I. For this reason we calculate the projection 

+1 

d(cose)0k(r)P,(cose) = 2i'AiRi{r), (17) 
where we have taken into account that the Legendre polynomials are orthogonal 
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/ d(cos e)Pi (cos e)Pi> (cos e) = - — -6' 
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We decompose the scattering amplitude, 
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"''(®) = ?i:E^'*'(2/ + i)PKcose), (19) 
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and we project the scattered wave function (fTIjl as 
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d(cos e) (e^'^'- + f{e)^)Pi (cos e) = 2 ^' (j, {kr) + Q ^) ■ (20) 



Next, we equate the asymptotic forms of the projections H17|l and ((2()|l . and thus, instead of Eq. ((15|1 . we obtain 

Ai sin(fcr - ttI/2 + 6i) = sin(A:r - ttI/2) + Cik- — . (21) 

r 

We compare the terms proportional to e"*'^'" and e*'^'', respectively, and find that Ai = e'*' and 

G = -^e-^^'/2p«' -1], (22) 

which, due to Eq. ((19|l . again provides the correct result Hlf)|l. 

Although the problem discussed here looks pure academic it was discovered in the course of concrete calculations. 
To simplify the calculation of a correlation function where the scattering wave function enters, we used the form 1(13(1 
with Ai = e"*' as is given in many books. We were interested in the complete sum of partial waves, and we used 
Eq. p3(l instead of Eq. 1(14(1 to exploit the orthogonality of Legendre polynomials. Needless to say the calculation 
went wrong, showing that the asymptotic expressions must be treated very carefully. 
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